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ABSTRACT 

We examine the properties of 18 long period (80-210 days) and very luminous (median absolute magnitude 
of M/ = -7.86 and My = -6.97) Cepheids to see if they can serve as an useful distance indicator. We find that 
these Ultra Long Period (ULP) Cepheids have a relatively shallow Period-Luminosity (PL) relation, so in fact 
they are more "standard candle"-like than classical Cepheids. In the reddening-free Wesenheit index, the slope 
of the ULP PL relation is consistent with zero. The scatter of our sample about the W/ PL relation is 0.23 mag, 
approaching that of classical Cepheids and Type la Supernovae. We expect this scatter to decrease as bigger 
and more uniform samples of ULP Cepheids are obtained. We also measure a non-zero period derivative for 
one ULP Cepheid (SMC HV829) and use the result to probe evolutionary models and mass loss of massive 
stars. ULP Cepheids main advantage over classical Cepheids is that they are more luminous, and as such show 
great potential as stellar distance indicators to galaxies up to 100 Mpc and beyond. 
Subject headings: Cepheids — stars: distances — stars: mass loss — distance scale 



1. INTRODUCTION 

A reliable method of measuring the physical distance to 
astrophysical objects has always been sought after in obser- 
vational astronomy (e.g., Bessel 1839). In the era of "pre- 
cision cosmology", the need for accura te physical distance 
measurements has been amplified (e.g ., ISpergel et al.ll2003l 
iRiess et all 12004 iTegmark et"aTJ 120041) . Accurate and pre- 
cise distance indicators hold the key to pinning down the 
value of the Hubble constant {Ho) and many other cosmo - 
logical parameters (see discussion in, e.g jMacri et al.l l2006). 
A number of methods have been employed to determine ex- 
tragalactic distances, with varying degree of success (e.g., 
Freedman et al. 2001). The construction and reliability of the 
"cosmological distance ladder" depends crucially on Cepheid 
variables being able to provide precise and accurate distances 
to nearby (d < 20 Mpc) galaxies. The quest for such distances 
has been an arduous journey for almost a hundred years, with 
ma ny dramatic twists and turns (for a review of early years, 
see iBaadal 19561 for a recent review see, e.g., Macri 2005). 

Cepheids offer several advantages as distance indicators. 
Massive stars (> 5M Q ) make an excursion through the in- 
stability strip and most, if not all, of them become Cepheid 
variables. These variable stars are relatively bright (My ~ -4 
for a P ~ 10 day Cepheid) and often have large brightness 
variations (amplitude ~ 1 mag) with a characteristic "saw- 
tooth" shaped light curve. Their intrinsic brightness, com- 
bined with their light curve shape and colors, make them easy 
to distinguish from other classes of variable stars. As a re- 
sult, Cepheids have been detected and studied in a signifi- 
cant number of star-forming galaxies. The physical mecha- 
nisms underlying Cepheid pulsation are well understood, in- 
cluding th e observed tight period-luminosity (PL) relation- 
ship (e.g., Chiosi et al. 1993). The small scatter in the PL 
relation allows d istance measurements precise to ~ 5% (e.g., 
iMacri et aL 2006). For these reasons, Cepheids are commonly 
used to calibrate other distance indicators, forming the base of 
the cosmological distance ladder. 

Despite their many advantages as a distance indicator, 
Cepheid distances also have some shortcomings. Most 



Cepheids have an intrinsic brightness of My > -5, so with 
the current instrumentation they can be only used to mea- 
su re distances to < 30 M pc (the largest Cepheid distance 
in iFreedman et al.1 [2001 is ~ 22 Mpc). Observations of 
Cep heids in distant galax ies are also hindered by blend- 
ing (Mochejskaet al. 2000) — as young stars, Cepheids live 
in close proximity to the crowded star-forming regions of 
their host galaxies, and are thus likely to have another star 
of similar brightness on the scale of a typical instrumen- 
tal point-spread-function (PSF). The effect of blending be- 
comes worse as the squa re of the distance to the host galaxy 
(Stanek & Udalski 1999), again limiting the usefulness of 
Cepheids to measuring distances < 30 Mpc even with high 
resolution instruments such as the Hubble Space Telescopes 
(HST). Ideally, we would like to find a distance indicator that 
shares the good properties of classical Cepheids, but is even 
more luminous, allowing us to observe it further away and be 
less susceptible to blending. In this paper we discuss such a 
possible distance indicator, namely Ultra Long Period (ULP) 
Cepheids. 

We define ULP Cepheids as fundamental mode Cepheids 
with pulsation periods greater than 80 days. Several such 
Cepheid s have been a lready discovered in the pioneering 
study of lLeavittl (119081) . However, ULP Cepheids have tra- 
ditionally been ignored for distance measurements as they 
are PL outliers. Indeed, the observed PL relation flat- 
tens for Cepheids with periods greater t han 100 days (e.g. , 
iGrieve et alJl985tlFreedrnan et alJ1992h . lGrieve et al] ( 1 1985b 
suggests that long period Leavitt Variables could be used for 
distance measures — unfortunately that idea has not perme- 
ated through the community. We argue that the flattening of 
the PL at long periods actually improves the usefulness of 
ULP Cepheids as distance indicators because it makes them 
a good standard candle in the traditional sense. We note sev- 
eral additional advantages of ULP Cepheids over lower period 
Cepheids due to their increased luminosity. ULP Cepheids 
could be used as a stellar distance measure to the Hubble Flow 
(up to ~ 150 Mpc) — several times the current observational 
limit. In Section [2] we describe our sample compiled from 
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the literature. The ULP Cepheid PL relation is discussed in 
Section |3] Section H demonstrates how ULP Cepheids may 
provide the additional benefit of testing massive stellar evolu- 
tionary models. We summarize our results in Section[5] 

2. SAMPLE 

We have assembled a sample of ULP Cepheids from the lit- 
erature and list their reported positions, periods, and mean V 
and / magnitudes (see TableQ}. We adopt the periods, redden- 
ing values, and distance moduli found in these sources except 
in the case of the Magellanic Clouds (see below). Our pri- 
mary criteria for selecting the sample was the existence of V 
and / data calibrated on the standard Johnson/Kron-Cousins 
magnitude system using Landolt standards (with the possible 
exception of the Magellanic Clouds; see below). The ULP 
distinction is applied to fundamental mode Cepheids with pe- 
riods greater 80 days. We combed the recent literature for 
reports of such variable stars. 

Magellanic Clouds Our sample includes four LMC and 
three SMC ULP Cepheids. IFreedman et all (Il985h calibrated 
photoelectric observations of these Cepheids and transformed 
them to the Johnson/Kron-Cousins standard system. The 
m ean flux-weighte d photo metry for the six Cepheids reported 
in IFreedman et al.l (1 19851) agrees with the L andolt standard 
star calibrated results of uvloffett et ail d!998t) to within 0.04 
mag, suggesting that the standard photometric system calibra- 
tion is robust. Mean flux- weighted photometry for HV1956 
is obtained from th e All Sky Automa ted Survey (AS AS; 
iPoimanskil 19971) and lMoffett et all dl998l) . The V light curves 
of six of these ULP Cepheids were obtained from AS AS. 
HV2883 was not targeted by AS AS, and it s V light curve 

J hotom etry was obtained fromlMadord (.1 975), van Genderen 
1981 . and iMoffett et alj dI998j). jMoffett et al. ( 1998} pro- 
vide the / light curve data for the ent ire sample. We applied 
the a nalysis of variance technique (Schwarzenberg-Czerny 
1989) to the seven Harvard Variable light curves in FigureQ]to 
obtain the periods listed in Table Q] We adopt total reddening 
values of E(B-V) = 0.14 mag a nd E(B-V) = 0.09 mag for 
the LMC and SMC, respectively dUdalskiet al.ll 19991) . We as- 
sume a distance modulus (DM) of (m— M)q =18.5 mag to the 
LMC for consistency with the sources li sted in Ta ble Q] The 
SMC DM used is (m- M) =18.93 mag dHilditch et alT2005l; 
iKeller & Woodl2006l) . The LMC (SMC) hosts ULP Cepheids 
with periods of 98.6, 109.2, 118.7, and 133.6 (84.4, 127.5, 
and 210.4) days. The LMC has gas phase oxygen abundance 
12+log(0/H) = 8.39±0. 10 (Pagel et al. 197 8) while the SMC 
is 12 + log(Q/H) = 7.98 ± 0.10"(Peimb ert & Torres-Peimbertl 
fl976h . 

The Araucaria Project (Pietrzvris ki et al.l 12002) is a photo- 
metric survey of Local and Sculptor Group galaxies and their 
Cepheid populations. The primary goal is to more accurately 
determine the distances to these galaxies and to character- 
ize the dependence of various stellar distance indicators on 
metallicity and age. The Araucaria Project has observed ULP 
Cepheids in the following galaxies. 

NGC 55 Five ULP Cepheids were found in NGC 55 
(Pietrzvrisk i et alj|200 6). Observations were taken with the 
Optical Gravitation Lensing Experiment (OGLE) detector on 
the Warsaw 1.3 m telescope at Las Campanas Observatory. 
They estimate that the calibration procedure used to transform 
their photometric data from the OGLE filters to the standard 
system produced errors < 0.03 mag. Follow up observations 
in the IR revealed a total reddening of E(B-V) = 0.13 mag 
dGieren et alj|2008l) . Their multi-wavelength PL analysis pro- 



duced a DM to NGC 55 of 26.43 ± 0.04 ± 0.08 mag (statisti- 
cal and systematic errors, respectively). NGC 55 hosts ULP 
Cepheids with periods of 85.1, 97.7, 112.7, 152.1, and 175.9 
days. The oxygen a bundanc e of NGC 55 is 12 + log(0/H) = 
8.05 ±0.10 (Tu llmann et ai] l2003h. 

NGC 6822 One ULP Cepheid was found in NGC 6822 
dPietrzvriski et al.N2004l) . The filters and telescope u sed are 
identical to those of NGC 55 dPietrzvriski et al. 2006). Sim- 
ilarly, the reported calibration error onto the standard system 
is < 0.03 mag. As in the multi-w avelength follow up study 
of NGC 6822 ( Gierenetai1l2006l) . we adopt a total redden- 
ing of E(B-V) = 0.36 mag. The lone ULP Cepheid i n NGC 
6822 has a period of 123.9 days. iGieren et alJ (12006) calcu- 
late a DM to NGC 6822 of 23.31 ±0.02 ±0.06 mag (statistical 
and systematic errors, respectively). NGC 6822 has a similar 
oxygen abundance to NGC 55 of 12 + log(0/H) = 8.11 ±0.10 
(Pei mbertetaT1l2005l). 

NGC 300 IGieren et all ( 120041) found three ULP Cepheids 
in NGC 300. Again, OGLE filters were used for the obser- 
vations. Their calibration onto the standard system has a re- 
port ed error < 0.03 mag . A multi-wavelength study of NGC 
300 dGieren et al.ll2005b determined a reddening-free DM of 
26.37 ±0.05 ±0.03 mag (statistical and systematic, respec- 
tively) using a total reddening of E(B-V) = 0.10 mag. ULP 
Cepheids of 83.0, 89. 1, and 1 15.8 days are observed in NGC 
300. NGC 300 has a strong metallicity gradient; therefore we 
adopt mean Cepheid radial distance of 4 kpc and apply the 
averaged gradient of Urbanej a et alj d2005l) to obtain a mean 
oxygen ab undance value of 1 2 + log(Q/H) = 8.25 ± 0.22. 

I Zw 18 AloisietalJ d2007l) discovered three ULP Cepheids 
from the extremely metal poor galaxy I Zw 18, though 
they could not confirm one candidate. A follow-up study 
dFiorentino et al . 2008) presents flux weighted mean photom- 
etry but no data, so the light curves of these objects could 
not be included in Figure Q] The ULP Cepheids have peri- 
ods of 129.8 and 125.0 days (Tabl ed). After accounting for 
extinction E(B-V) = 0.032 mag. lAloisi et all d2007l) use the 
red gi ant branch tip to determ ine a DM of 31.30 ±0.17 mag 
while lFiorentinoetal.1 d2008l) find a DM of 3 1 .35 ± 0.26 mag 
via pulsation models. We use the former measurement as it 
is considered more reliable by the authors. We do not include 
these two Cepheids in the upcoming PL determination as I Zw 
1 8 is a full dex more metal poor than th e other galaxies in this 
samp le (12 + log(0/H) = 7.2 ± 0.10; ISkillman & Kennicuttl 
1993). There is an increasing amount of support for a metal- 
licity dependent PL (e.g. Sandage et al. 2008) and including 
these Cepheids in our PL analysis would greatly increase the 
metallicity dispersion of the host galaxies in our sample. We 
do, however, make use of them to examine the ULP PL rela- 
tion dependence on metallicity. 

2.1. Absolute Photometry 

The ULP Cepheid sample and its mean, flux-weighted pho- 
tometry in the standard system can be found in Table Q] We 
assume that the photometric error associated with each ULP 
Cepheid is negligible when compared to the intrinsic scatter 
of the PL relation. We transform these measurements to ab- 
solute magnitudes via: 



Mi = mi-DM-A h i = I,V 



(1) 



where M, is the absolute magnitude in either the V or /, w,- is 
the apparent magnitude; DM is the reddening free distance 
modulus; and A, is the extinction in the V or /. We use 
the extinction law A v = 3. 24E(B-V) andA 7 = l.96E(B-V) 
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FIG. 1 . — The V (open circles) and / (filled circles, where available) band light curves of the ULP Cepheids. Each panel spans 2.4 magnitudes along the y-axis. 
The phase is given along the x-axis. The Cepheid identification (see column 1, Table[T} is listed in the upper left of each plot while the period (in days) is in the 
upper right. 



dSchlegel et al.|[l998h . We d efine the Wesenheit magnitudes 
as: W, = I-1.55(V-l) (e.g. JUdalski et al1fl999h . 

The color-magnitude diagram (CMD) highlights several 
important characteristics of the ULP Cepheid data set (Fig- 
ure |2). ULP Cepheids are the luminous counterparts to 
shorter period Cepheids in color magnitude space. Our sam- 
ple clearly populates the luminous region of the instabil- 
ity strip. Future Cepheid studies can use ULP Cepheids to 
push Cepheid distance measurements well beyond the current 
~ 30 Mpc limit as our sample's median absolute magnitude is 
M/(My) = -7.86(-6.97) (see SectionEJ. The intrinsic bright- 
ness of ULP Cepheids makes them ideal candidates for dis- 
tance indicators to galaxies where classical Cepheids cannot 
be observed. 

3. DISTANCE MEASUREMENTS WITH ULP CEPHEIDS 

The Cepheids in our sample have been ignored as dis- 
tance indicators as they do not follow to the standar d period- 
luminosity relationship (e.g., Freedma n et al.lll992b . In this 
section, we examine the characteristics of our ULP Cepheid 



sample in the period-luminosity plane and explore their via- 
bility as a distance indicator. We determine PL relations of 
our sample in the V, I, and W/ in Section [3~T| while the metal- 
licity dependence of our results is presented in Section [3~2l 

3.1. Period Luminosity Relations 

Using the data in Table [1] we construct PL diagrams in V, 
I, and W/ (Figures [3] IU [5] respectively). In each case, ULP 
Cepheids are compared to fundamental mode SMC Cepheids 
(hereafter, this control sample will be referred to as OGLE 
Cepheids). The OGLE Cepheid sample contains over 1100 
fundamental mode Cepheids with periods ranging from 0.5 
to ~ 50 days; however, we only plot the 70 Cepheids with 
periods of greater than 10 days. To quantify our compari- 
son, we perform a linear fit on both samples in each PL dia- 
gram: the slopes of the OGLE Cepheid PL relations (here- 
after PL^mc; dotted lines in the figures) are adopted from 
Udals ki et al.l (Q999) while the intercepts are chosen to mini- 
mize chi-square. We employ linear least square fitting of the 
ULP Cepheid sample to identify their PL relation (hereafter 
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TABLE 1 
Ultra Long Period Cepheids 



ID 


Host Galaxy 


RA 


DEC 


P 


<V> 


(V-I) 


w, 


E(B-V) 


(m-M)o 


Vo 


(V-/)o 


W,o 


12 + log(0/H) 


References 






(J2000.0) 


(J2000.0) 


(days) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


dex 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 



LMC HV883 
LMC HV2447 
LMC HV2883 
LMC HV5497 
SMC HV1956 
SMC HV821 
SMC HV829 
NGC 6822-1 
NGC 55-1 
NGC 55-2 
NGC 55-3 
NGC 55-4 
NGC 55-5 
NGC 300-1 
NGC 300-2 
NGC 300-3 
IZw 18-1 
IZw 18-2 



LMC 
LMC 
LMC 
LMC 
SMC 
SMC 
SMC 

NGC 6822 
NGC 55 
NGC 55 
NGC 55 
NGC 55 
NGC 55 
NGC 300 
NGC 300 
NGC 300 
IZw 18 
IZw 18 



05:00 
05: 19 
04:56 
04:55 
01 : 04 
00:41 
00:50 
19:45 
00 : 14 
00 : 15 
00 : 14 
00 : 15 
00 : 15 
00:55 
00:54 
00:54 



08.6 
31.4 
27.0 
40.0 
15.9 
43.5 
28.9 
02.0 
13.0 
12.0 
36.6 
14.3 
10.1 
11.6 
35.0 
54.3 



-68 : 27 : 03 
-68:41 : 12 
-64 : 41 : 26 
-66 : 25 : 39 
-72 : 45 : 20 
-73 : 43 : 24 
-72 : 45 : 09 
-14 : 47 : 33 
-39 : 08 : 42 
-39: 12: 18 
-39 : 1 1 : 09 
-39: 13: 17 
-39 : 12 : 26 
-37 : 33 : 55 
-37 : 35 : 01 
-37 : 37 : 02 



133.6 
118.7 
109.2 
98.6 
210.4 
127.5 
84.4 
123.9 
175.9 
152.1 
112.7 
97.7 
85.1 
115.8 
89.1 
83.0 
129.8 
125.0 



12.12 
11.99 
12.41 
11.92 
12.28 
11.92 
11.97 
17.86 
19.25 
19.56 
20.18 
20.54 
20.84 
20.13 
19.71 
19.26 
23.56 
23.47 



1 .09 

1.12 
1.07 
1.11 
0.83 
1 .03 
0.91 
1.40 
0.84 
0.95 
1.05 
1.25 
1.38 
0.97 
1.02 
0.77 
0.74 
0.91 



9.34 
9.13 
9.68 
9.09 
9.95 
9.29 
9.65 
14.29 
17.11 
17.14 
17.51 
17.35 
17.32 
17.66 
17.12 
17.30 
21.67 
21.15 



0.14 
0.14 
0.14 
0.14 
0.09 
0.09 
0.09 
0.36 
0.13 
0.13 
0.13 
0.13 
0.13 
0.10 
0.10 
0.10 
0.03 
0.03 



18.50 
18.50 
18.50 
18.50 
18.93 
18.93 
18.93 
23.31 
26.43 
26.43 
26.43 
26.43 
26.43 
26.37 
26.37 
26.37 
31.30 
31.30 



-6.83 
-6.96 
-6.54 
-7.03 
-6.94 
-7.30 
-7.25 
-6.60 
-7.60 
-7.28 
-6.67 
-6.31 
-6.01 
-6.55 
-6.97 
-7.42 
-7.84 
-7.93 



0.91 

0.94 
0.89 
0.93 
0.71 
0.92 
0.80 
0.94 
0.68 
0.79 
0.88 
1.08 
1.22 
0.85 
0.89 
0.65 
0.70 
0.87 



-9.16 
-9.37 
-8.82 
-9.41 
-8.98 
-9.64 
-9.28 
-9.02 
-9.33 
-9.29 
-8.92 
-9.08 
-9.12 
-8.71 
-9.25 
-9.07 
-9.63 
-10.15 



8.39±0.12 
8.39±0.12 
8.39±0.12 
8.39±0.12 
7.98±0.10 
7.98±0.10 
7.98±0.10 
8.11 ±0.10 
8.05 ±0.10 
8.05 ±0.10 
8.05 ±0.10 
8.05 ±0.10 
8.05 ±0.10 
8.25 ±0.22 
8.25 ±0.22 
8.25 ±0.22 
7.21 ±0.10 
7.21 ±0.10 



1,2,3 

1,2,3 

1,2,3 

1,2,3 

5,2,3;4 

1,6,3;4 

1,6,3;4 

7,8,9 

10,11,12 

10,11,12 

10,11,12 

10,11,12 

10,11,12 

13,14,15 

13,14,15 

13,14,15 

16,17 

16,17 



REFERENCES. — References: l. | Freedman et alj i 1985); 2. Pagel et al. 1 1978); 3. Udalski et al. 1 1999); 4. Hilditch et al. i2005) and Keller & Wood .2006); 5. Ppjmanski . 199 7], AS AS survey) ; 
6. Peimbert & Torres-Peimbert 1 1976); 7. Pietrzyriski et al. i 2004); 8. Peimbert et al. (2005); 9. Gieren et al. i 2006) 10. Pietrzyriski et al. (2006); 1 1. Tullmann et al. i 2003); 12. Gieren et al. |2008) 
13. Gieren et al. i 2004); 14. Urbaneja et al. (2005); 15. Gieren et al. 1 2005); 16. Fiorentino et al. 1 2008); 17. Skillman & Kennicutt 1 1993). 

NOTE. — ULP Cepheids in our sample grouped by host galaxy. (1): Cepheid Identification. (2): Host galaxy name. (3,4): Right Ascension and Declination in J2000 coordinates. (5): Period 
in days. (6): Apparent mean V magnitude of Cepheid. (7): Apparent (V — I) color of Cepheid. (8): Apparent Wesenheit magnitude (definition in text). (9): Reddening towards host galaxy. (10): 
Reddening-free distance modulus. (11): Absolute V magnitude. (12): Absolute (V —I) color. (13): Absolute Wesenheit magnitude. (14): Metallicity: 12 + log(0/H). (15): First reference is for 
photometry, reddening, and distance modulus (except where noted in text). Second reference refers to metallicity. If third reference is present, its reddening and distance modulus measurements 
supercedes the first. 



PLulp)- We omit errors in distance moduli and extinction in 
this demonstration as they are small when compared to the 
overall scatter of the sample. The parameters of these fits and 
the RMS of each data set in V , I, and W/ are listed in Table [2] 
Despite our sample ranging in period from 83 to 210 days, 
ULP Cepheids occupy a small region of luminosity space. We 
note that the ULP Cepheids from I Zw 1 8 are not included in 
this analysis for reasons outlined in Section[2] 

The V PL diagram is Figure [3] The PL,smc fit has a slope 
of -2.76 and the RMS of the OGLE Cepheid sample about 
this fit is 0.25 mag. The ULP Cepheid sample has 76% more 
scatter (RMS=0.44 mag) about the PL SMC fit. This discrep- 
ancy in scatter has led to the standard pract ice of removing 
ULP Cepheids from PL relation studies (e.g., Freedm an et alj 
1985) and the significant increase in RMS suggests that the 
ULP Cepheids do not conform to the classical PL relation. 
If we determine the PL relation of ULP Cepheids alone we 
find PL ULP has a slope (-1.09 ±0.94) that is flatter than 
PLsmc (though the slopes are within 2er of each other, see 
Tabled The RMS of our sample to PLulp is 0.40 mag; only 
marginally better than the ULP Cepheid scatter about the es- 
tablished PLsmc relation. In V, the ULP Cepheid sample does 
not follow a statistically distinct and unique PL. 

The longer the wavelength the less reddening is a con- 
cern. The accuracy of distance measurement with Cepheids 
increases going from V to I (PL diagram in Figure|4]i. PL^mc 
has a slope of -2.96 and the OGLE Cepheid sample's RMS is 
0.19. The ULP Cepheid scatter about this fit is 116% larger 
(0.41 mag). The RMS of the ULP Cepheids is reduced to 0.31 
mag when using the PLulp fit (slope= -0.57 ±0.73). PLulp is 
approximately five times as flat as PLsmc and the two slopes 
are distinct at the 3a level. While ULP Cepheids show the 
same general trends with regards to period, luminosity, and 
color as normal Cepheids, significant statistical differences 
between the two populations are apparent in the /-band PL 



relation. 

To further reduce the uncertainty associated with red- 
dening in our PL analysis, we repeat the procedure us- 
ing the reddening-free Wesenheit Index (W/) introduced in 
Mado re& Freedmanl d!991h . The W/ PL diagram illuminates 
the advantages of this reddening-free index for Cepheid dis- 
tance measurements (Figure |5}. The OGLE Cepheid sample 
has a very tight relation between period and luminosity; with 
small scatter about the PLsmc fit (slope of -3. 28; RMS of only 
0.12 mag). The ULP Cepheid RMS about PL W c is 0.47 mag. 
The 4-fold increase in scatter implies the ULP Cepheid and 
the OGLE Cepheid samples do not conform to the same PL 
relation. The PLulp fit slope is flatter than its PL^mc coun- 
terpart at the 6a level (-0.05 ± 0.54 vs. -3.28) and the ULP 
fit slope is consistent with zero slope. The scatter of the ULP 
Cepheids is only 0.23 mag about the PL^p relation. This 
scatter is still 92% more than that of the OGLE Cepheid sam- 
ple; however, the ULP sample is relatively small and hetero- 
geneous. We note that the scatter of ULP Cepheids about the 
PLulp fit is smaller than that of the OGLE Cepheid sample 
about the nominal PL relation in V and on par with the same 
in/. 

Several trends in PL space are apparent as one examines 
the ULP Cepheid sample in V, then /, and finally in the W/ 
index. As reddening is reduced, the PLulp parameters are in- 
creasingly different from those of PLsmc- The PLulp fit is 
more shallow as one moves from V to W/. In the W/, the 
PLulp relation reveals that ULP Cepheid luminosity becomes 
statistically independent of period. In essence, ULP Cepheids 
behave as bright, standard candles in the reddening-free in- 
dex. 

3.2. Metallicity Dependence 

An uncertainty of the Cepheid PL and its derived distance 
measurement is its sensitivity to the metallicity of the stars 
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(e.g.. lFreedman&Madodll990HKennicuttetalJ ll998). Our 
sample of ULP Cepheids contains six different host galax- 
ies that span a range of ~ 1.2 dex in 12 + log(0/H) from 
7.22 to 8.39, providing an opportunity to investigate the de- 
pendence of the ULP Cepheid PL on metallicity. We plot 
the residual of each ULP Cepheid to the PLulp.w, fit listed 
in Table [2] as a function of metallicity (Figure |6). In other 
studies, linear fits of PL residuals have determined a correc- 
tion fa ctor, 7, between 0.0 and -0.4 mag dex" 1 (see Figure 
1 of iRomaniello et alj 120081) . Recently. iMacrietal J (12006b 
used the metallicity gradient in NGC 4258 to determine 7 = 
-0.29 ±0.09 ±0.05 mag dex" 1 (random and systematic errors, 
respectively). We overlay this relation (normalized to our data 
set) in Figure[6]for reference. If the I Zw 18 Cepheids are con- 
firmed, it suggests a stronger correlation between PL offset 
and metallicity than is evident in lower period Cepheids (e.g., 
iKochanekl [l997t: iKennicuttet alj|1998l) . However, we note 
that we do not take into account any reddening or DM errors 
in this analysis. As such, we do not claim a specific metal- 
licity dependence for ULP Cepheids. We simply demonstrate 



TABLE 2 
Least Square Fit Values: y - b + a * 
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PL SMC.V 
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PL SMC.I 
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-0.57 ± 0.73 
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Period-Luminosity: Wj 
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Period-Luminosity: Wj 


index 


ULPs 
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-9.06 ± 1.12 


-0.05 ± 0.54 
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Note. — The fit values oflhc PL relationships in V. I. and W/. For cadi photometric system, > 
of classical SMC (PL SMC ) and ULP (PL(/ LP ) Cepheids. 



9 Ihc PL relation 



that the residuals to the PLuip fit are broadly consistent with 
the range of values presented in the literature to date. 

At this time we do not apply a metallicity correction to our 
ULP Cepheid PL relations. Precise gas phase oxygen abun- 
dance measure ments are difficult to obtain (for a review see 
Bresolin 2006) and we adopt a minimum metallicity error of 
0.1 dex. The ULP Cepheid sample must grow in size and the 
PL analysis must be more detailed to determine if I Zw 18 
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FIG. 3. — The V Period Luminosity relationship for the OGLE SMC Cepheids (black dots) and ULP Cepheids (open symbols). The dashed line is the PL 
relation adopted from Udalski et al. ( 1999), with a slope of —2.76. The least square fit of the ULP Cepheid subsample yields a flatter slope of —1.09 ±0.94 (red 
line) and the RMS is 0.40 mag. The residuals to the PLsmc fit are shown in the bottom panel (black, open squares). Residuals to the PL(/u> are given for the ULP 
Cepheid sample (red symbols). 

is an exception and to characterize the functional form of the 
ULP Cepheid PL sensitivity to metallicity. 



4. USING ULP CEPHEIDS TO PROBE THE EVOLUTIONARY 
MODELS OF MASSIVE STARS 

Most Cepheid variables cross the instability strip three 
times (e.g., Bono et al. 2000; Pietrukowicz 2001). One can 
determine which crossing a Cepheid is undergoing by mea- 
suring its period change, dP/dt. The first and third crossings 
are associated with positive dP/dt while Cepheids exhibit a 
decreasing period on their second crossing. We investigated 
the light curves of our sample for signs of period changes 
by comparing photometry taken over the last 30 years. One 
ULP Cepheid, HV829, exhibits a negative period change of 
about 1.5 days. For this Cepheid, we compiled pho tomet- 
ric data taken during 1970-1976 from iMadorel d 1975b and 
Ivan Gender en ( 1983) and 2000-2004 data from the ASAS cat- 
alog (See Figure |7]). This result is confirmed by a measured 
f eriod of 87.63 days i n iPavn e-Gaposchki n & Gaposchkir] 
1966) and 85.2 days by MoffettetalJ il998j), firmly estab- 



lishing HV829 as a Cepheid on its second crossing. No other 
Cepheid in our sample exhibited a measurable period change. 
ULP Cepheids occupy a mass range that is ideal to probe 
high mass evolutionary models along the instability strip. We 
plot the CMD of our ULP Cepheid samp l e and overlay the 
evolutionary tracks of lLejeune & SchaererJ (12001) (Figure [8]). 
The location of ULP Cepheids in the color magnitude dia- 
gram suggest masses between 13 and 20M Q , depending on 
the assumed metallicity (here, we choose SMC metallicity). 
Our sample clearly probes a mass range unexplored by cur- 
rent Cepheid studies. Note the evolutionary tracks at 15 and 
2OM only cross the instability strip once, regardless of as- 
sumed mass-loss rate. However, we have shown and con- 
firmed that HV829 is undergoing its second crossing; a result 
at odds with the models. We note that the evolutionary model 
represents some "mean" behavior based on assumptions of 
stellar chemical composition, overshooting, and other param- 
eters. Nevertheless, our result is one of the few observational 
constraints on high mass stellar evolutionary models. HV829 
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FIG. 4. — The / Period Luminosity relationship for the OGLE SMC Cepheids (black dots) and ULP Cepheids (open symbols). The dashed line has a slope of 
—2.96 and is the PL relation adopted from Udalski et al. ( 1999). The least square fit to the ULP Cepheid subsample produces a flatter slope of —0.57 ± 0.73 (red 
line) and the RMS is 0.3 1 mag. The residuals to the PLsmc fit are shown in the bottom panel (black, open squares). Residuals to the PL(//j> are given for the ULP 
Cepheid sample (red symbols). 

may indicate that massive stars behave differently than ex- 
pected. Future stellar evolutionary models in this mass regime 
should take this into account. 



5. DISCUSSION 

We have presented, for the first time, a collection of Ul- 
tra Long Period (P > 80 days) Cepheids from the literature 
and demonstrated their viability as a distance indicator. In the 
past, ULP Cepheids have been ignored as distance indicators, 
and in fact many stellar variability searches did not extend 
their cadence and search strategy to allow for their discovery. 

In V, I, and especially reddening-free W/ magnitude, ULP 
Cepheids have a relatively flat PL relatio n compared to the re - 
spective relations for classical Cepheids (Udalski et al. 1999). 
The dispersion in both the classical and ULP Cepheid popula- 
tions about their respective PL relationships becomes smaller 
as one moves from V to / to W/, and the discrepancy be- 
tween the slopes increases. Most notably, the slope of the 
ULP Cepheid W/ PL relation is significantly shallower than 



the standard SMC PL relation (slope is -0.05 vs. -3.28). The 
reddening-free Wesenheit index produces the tightest ULP 
PL relation, with RMS r esidual of only 0.2 3 mag (See Fig- 
ure|5]l. Other papers (e.g. Kan bur et al.l2 007) have found non- 
linearites in the PL relationship at lower periods. However, 
the change in slope seen here is much more dramatic and it 
is unlikely that it shares a physical connection with the PL 
changes at lower periods. 

Our ULP PL scatter in W/ is already less than that of the 
initial peak brightness vs. absolute magnitude relation for 
Type la Supernova (~ 0.3 mag; Phillips (1993)). A huge 
amount of effort has been necessary to increase the number of 
observed Type la SNe, refine the calibration technique (e.g., 
iPrieto et al1l2006f) . and to obtain t he low 015-0.2 mag scat- 
ter in the relation seen today (e.g jjha et al.11200 7). We expect 
future observational and theoretical studies of ULP Cepheids 
will further decrease the already fairly small scatter found in 
this first analysis. 

We strove to find all the known ULP Cepheids in the liter- 
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FIG. 5. — The absolute Wesenheit magnitude Period Luminosity relationship for the OGLE SMC Cepheids (black dots) and ULP Cepheids (open symbols). 
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The residuals to the PLj M c fit are shown in the bottom panel (black, open squares). Residuals to the PL(//.p are given for the ULP Cepheid sample (red symbols). 

a period of ~ 121 days) at 100 Mpc and the brightest ULP 
Cepheids at ~ 150 Mpc in a reasonable amount of time. Since 
the luminosity of an ULP Cepheid is a weak function of its pe- 
riod, relatively accurate distances would not require as precise 
period measurements as is needed for classical Cepheids. We 
encourage future variability proposals to search for Cepheids 
with periods greater than 100 days. 

We note two concerns in using ULP Cepheids for distance 
measurements. It is obvious from our sample size that ULP 
Cepheids are far less common than classical Cepheids. The 
relatively small ULP Cepheid population will make preci- 
sion distance measurements less practical. As the sample size 
grows and the ULP Cepheid PL relations become established, 
a single ULP Cepheid observation may yield distance mea- 
surements accurate to 10-20%. Blending is always a con- 
cern in Cepheid studies. Blending can compromise Classical 
Cepheid measurements as stars o f comparable bri ghtness are 
likely to lie within a single PSF (Mochejska et al. 2000), es- 
pecially at larger distances. We expect blending to be less of 



ature, but our sample is likely not a complete census of ULP 
Cepheids. We note that the Araucaria Project found ULP 
Cepheids in three of the five galaxies they observed at the 
time of this analysis, so the ULP Cepheid sample should grow 
at a reasonable rate as Cepheid studies are extended to more 
galaxies and previous data sets are reanalyzed with longer pe- 
riod searches. There is also evidence ULP Cepheids exist in 
a broad range of metallicities, as preliminary data analysis of 
a M8 1 variability survey with the Large Binocular Telescope 
has already discovered at least one ULP Cepheid (Kochanek, 
private communication). 

In addition to following a fairly tight PL relation, the ULP 
Cepheids are also very luminous, with a median absolute 
magnitude in I(V) for our sample of — 7.86(— 6.97). Us- 
ing WFPC3, HST can obtain 10% photometry at V = 28 or 
/ = 27 (DM = 35 — distance of 100 Mpc for the median ULP 
Cepheid) in about 10 orbits, while only two orbits are needed 
to reach a DM of 34. As one would only need a few orbits 
per epoch, one could detect the median ULP Cepheid (with 
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a problem for our ULP Cepheids even at large distances sim- 
ply because they are so bright and there are very few stars of 
comparable brightness in a given galaxy. The effect of blend- 
ing on ULP Cepheid observations will need to be investigated 
further as the sample is increased. 

The data set and analysis herein provides meaningful con- 
straints on theoretical work on Cepheids in this mass and pe- 
riod range. We have shown that ULP Cepheids show strong 
evidence for a different, flatter PL relation than their lower 
period cousins. Several papers have modeled Cepheid pul- 
sations and map ped these to theoretical PL relations (e.g., 
iBono et al1l2002l) . However, this work has not been reliably 
extended to the period range of our sample. Pulsation models 
in this period range would also help determine the intrinsic 
scatter to the PLj/l/ 5 relation. 

Our current sample is not large enough to constrain the sen- 
sitivity of the ULP Cepheid PL to metallicity. The PL resid- 
uals as a function of metallicity are consistent with the re- 
sults f or shorter per iod Cepheids (7 = -0.29 ±0.09 ±0.05 mag 
dex~': lMacri et al.ll2006l ). Note that the median 12+log(0/H) 



value for a galaxy in our sample is about 0.5 dex lower than 
the corresponding value for the HST Key Project Cepheid 
hosts. 

Period changes of ULP Cepheids have powerful implica- 
tions for stellar evolutionary models in this mass regime. 
We examined photometry from two epochs separated by 
thirty years to check for evidence of period changes in seven 
Magellanic Cloud ULP Cepheids. Only one, SMC HV829, 
showed a significant period change — from 85.9 ±0.3 to 
84.4 ± 0.4 days. The negative period derivative indicates that 
the Cepheid is undergoing its second crossing. Assuming the 
metallicity of the SMC, many current evolutionary models do 
not predict second crossings for Cepheids in this mass range 
(See Figure [§). Future models should incorporate this obser- 
vational result as it should place limits on several key input 
parameters. In addition to their potential as distance indica- 
tors, ULP Cepheids provide observational constraints on high 
mass stellar evolution models. 

To summarize, ULP Cepheids, while often dismissed in the 
past, are potentially a powerful distance indicator, probe of 
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phased with the period shown at the top of each pane. The pulsation period has clearly decreased, from 85.9 days to 84.4 days over approximately 30 years. 
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PL metallicity sensitivity, and also a probe of massive star 
evolutionary models. ULP Cepheids could provide the first 
direct stellar distance measurements to galaxies in the 50-150 
Mpc range, extending the cosmological distance ladder well 
into the Hubble flow. 
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